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Despite having entered the era of “precision cosmology,” the formation of galaxies within the favored CDM cos¬ 
mological paradigm remains problematic. By relating our N-body/smooth particle hydro-dynamical simulation 
to an extensive range of Galactic and extragalactic observations, we shed light on the formation and evolution 
of the Milky Way and other late type galaxies. In light of recent observations of the stellar populations of 
extragalactic thick disks, we examine the proposal that the thick disk formed in a high redshift period charac¬ 
terized by gas rich merging. We show such a scenario to be consistent with color observations. We then follow 
the evolution of structural parameters of the subsequently formed disk galaxy, from redshift ~1 to the present. 
Consistent with observation, little evolution in the ratio of scale-height to scale-length is found in our simulated 
galaxy, despite its somewhat chaotic origins. The simulated galaxy in this report forms part of a larger study 
of a suite of galaxies, with which these issues are to be studied in detail. 


1. INTRODUCTION 

The age of precision cosmology is upon us, with 
concordance values of cosmological parameters plac¬ 
ing us in an accelerating Universe, composed primar¬ 
ily of dark energy (^73%), and dark matter (^23%), 
with a mere smattering of baryonic matter (^4%). 
The evolution of structure on large scales is driven 
by gravity, and such structures as clusters and fila¬ 
ments of galaxies are exceedingly well explained by 
the ACDM paradigm. Yet on the scale galaxies, prob¬ 
lems persist. In some sense this is not surprising, as 
on smaller scales, baryonic processes become impor¬ 
tant; gas dynamics, star formation, quasars. Active 
Galactic Nucleii, supernovae, shocks, and ram pres¬ 
sure are some of the processes that add complexity 
beyond gravity driven structure formation. Our work 
is aimed at unravelling the major processes which de¬ 
termine how galaxies form within the context of hier¬ 
archical structure build. 

The recently emerging conviction that thick disks 
are prevalent in disk galaxies, and their seemingly 
ubiquitous old ages Q, means that the formation of 
the thick disk, perhaps more than any other compo¬ 
nent, holds the key to unraveling the evolution of the 
Milky Way^and indeed all disk galaxies. In an ear¬ 
lier study [ 3 , we proposed that the thick disk was 
formed in an epoch of gas rich mergers, at high red¬ 
shift. This hypothesis was based on comparing N- 
body/SPH simulations to a variety of Galactic and 
extragalactic observations, including stellar kinemat¬ 
ics, ages and chemical properties. Recently, the stellar 
populations in four nearby edge-on disk galaxies were 
resolved , and it was confirmed that thick disks are 
common features of spiral galaxies. Thick disk popu¬ 
lations of all galaxies in this sample are old and rel¬ 
atively metal rich, and no significant color gradient 
with height above the plane {AV—I/AZ) is apparent. 
These observations will provide further constraints to 


our thick disk formation scenario, which we probe in 
this study by examining how the ages and metallic- 
ties of our simulated thick disk stars vary with height 
above the plane. 

Meanwhile, in a study of 34 edge-on disk galaxies 
in the Hubble Deep Field 0, selected for apparent di¬ 
ameter larger than 1.3” and unperturbed morphology, 
it was shown that distant and local disk galaxies are 
globally similar in their relative thickness or flatness, 
as characterised by the ratio of scale-height to scale- 
length {hz/h\). Disks at redshift z ^ 1 are smaller, 
in absolute value, than present day disk galaxies, and 
have a flatness ratio, h^jhx, only slightly larger, by a 
factor ^1.5. We examine the evolution of our simu¬ 
lated disk galaxy in the epoch which follows thick disk 
formation, in light of these observations, to determine 
whether a choatic thick disk formation epoch is con¬ 
sistent with a galaxy whose structural parameters do 
not evolve much since around redshift 1. 

The galaxy presented in this conference proceeding 
is the first of an ongoing study of a suite of simulated 
galaxies, which will probe deeper into the issues pre¬ 
sented here, and will be detailed in 0 & |^. 

2. CODE AND MODEL 

We analyse data from the simulated late type 
galaxy from a previous study 0. Our Galactic 
Ghemodynamical code, (GCD+), models the effects of 
gravity, gas dynamics, radiative cooling, and star for¬ 
mation in a self-consistent manner. Type la and 
Type II supernova feedback are included, and we trace 
the lifetimes of individual stars when monitoring the 
chemical enrichment history of our simulated galax¬ 
ies. Details of GCD+ can be found in 0. We assume 
that 10^^ ergs is fed back as thermal energy from each 
SNe. An adiabatic phase is applied to gas heated by 
SNe II explosions. 
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Figure 1: Density plots of the evolution of the distribution of stars (upper panels) and gas (lower), shown face on, of 
the simulated galaxy between redshift z — l.Q and z = 1.1.Axis are 100 kpc. This epoch is characterised by a series of 
merger events. 
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Figure 2: Out of the plane direction (W) velocity 
dispersion versus log stellar age (Gyrs), plotted as red 
squares. Observations of solar neighbourhood stars are 
also shown, taken from (triangles with error bars) and 
m (which selects for thin disk stars, “+” symbols). An 
abrupt increase in velocity dispersion at ~10 Gyrs is 
apparent in the observations, and a similar abrupt 
increase is apparent in the simulated galaxy at ~8.5 
Gyrs. 


We use the semi-cosmological version of GCD+, with 
an isolated sphere of dark matter and gas as initial 
condition. Small scale density fluctuations are super¬ 
imposed on the sphere, parameterized by erg, which 
seed local collapse and subsequent star formation. 
Solid-body rotation is imparted to the initial sphere 
to incorporate the effects of longer wavelength fluctu¬ 
ations. For the model described here, the relevant pa¬ 
rameters include the total mass (5 x 10 ^^ Mq), baryon 
fraction ((2^ = 0.1), star formation efficiency, c*=0.05, 
spin parameter (A = 0.0675), and erg = 0.5. We em¬ 


ployed 38911 dark matter and 38911 gas/star parti¬ 
cles. 


3. RESULTS 

Our final simulated galaxy has properties resem¬ 
bling the Milky Way. The galaxy has a dominant, ex¬ 
ponential (scale-length^4.5 kpc), rotating disk com¬ 
prising primarily young, metal rich stars. It also has 
a low mass, low metallicity, pressure supported stel¬ 
lar halo (for details see , which examines the same 
simulated disk galaxy). 

3.1. Thick Disk Formation 

The epoch of thick disk formation is depicted in 
Figure [2| which shows density plots of the evolution 
of the distribution of star particles (upper panels) and 
gas (lower), shown face on, of the simulated galaxy 
between redshift 2 ; = 1.6 and z = 1.1. Axis are 100 
kpc. The velocity dispersion-age plot (Figure 0 pro¬ 
vides evidence that our simulated galaxy has a thick 
disk component, in close agreement with the velocity 
dispersion-age relation for solar neighbourhood stars 
0, using data of [l3|. The abrupt increase in velocity 
dispersion corresponds to end of the epoch depicted 
in Figure (21 This epoch, is by far the most chaotic 
period in the galaxies evolution. At the beginning of 
this epoch, at least four proto-galaxies of significant 
mass exist. These building blocks are gas rich, with 
combined gas mass of ^2.4x10^^ Mq, compared with 
stellar mass of r\j 7.0x10® Mq. By r\j 9 Gyrs ago, a 
single central galaxy has emerged. 
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Figure 3: the relation of the metallicity ([Fe/H]) of our 
simulated thick disk stars with height (kpc) above the 
plane. Very little gradient is apparent. 
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Figure 4: log age versus height for the thick disk stars of 
our simulated galaxy. There is no gradient apparent. 


Our thick disk stars have metallicity [Fe/H]^-0.5 
out to 2 kpc from the plane (Figure IHj). The nature of 
thick disk formation ensures that our simulated thick 
disk stars are old, and that there is little variation in 
age with |Z| (Figure 2|). 


3.2. Evolution of structural parameters 

The evolution of our simulated galaxy, in late 
epochs after the disk has formed, is shown in Fig¬ 
ure [5| where we plot the /-band luminosity at three 
times, z = 0.72 (left panels), ^ = 0.46 (middle), and 
z = 0 (right), both edge-on (upper panels) and face-on 
(lower). The scale parameters of our disks can be mea¬ 
sured by examination of radial and vertical luminosity 
profiles, shown in the /-band in Figure|n|and Figure[71 
respectively. It is evident that these parameters have 
not evolved much since the time that the disk galaxy 
has settled down. The luminosity surface density pro¬ 
files in Figure ini&H are shown at redshifts 2 ; = 0.72 in 
red, ^ = 0.46 in yellow and z = 0 in blue. We derive 
scale-lengths and scale-heights by fitting exponentials 


to these plots. The relative thickness to flatness can 
be parameterised by the ratio of scale-height to scale- 
length, hz/h\. This value is only larger by a factor of 
~1.2 in our simulated galaxy at z = 0.72 compared 
with the present. 


3 . 3 . Merger Histories 

We make now a statistical study of merger histo¬ 
ries in a ACDM model with Qq = 0.27, Aq = 0.73, 
and i/n = 71 kms“^Mpc“^, using a P^M algorithm 
[ 1 ^ . In our initial run, we simulated a comoving 
cubic volume of size Tbox = 100 Mpc, using 128^ par¬ 
ticles. The mass per particle was 1.801 x IO^^Mq. We 
identified a region of size 30 Mpc x 35 Mpc x 28 Mpc 
in which several Milky-Way-like halos formed. We 
zoom-in by increasing the resolution of the code by 2 
in each directions in the region of interest and rerun¬ 
ning the simulation. This simulation has 2.5 million 
particles, with mass per particle in the region of in¬ 
terest of 2.251 X IO^Mq. This is sufficient to resolve 
the relatively large building blocks which we are inter¬ 
ested in. We have 58 halos with masses in the range 
(0.5 — 1.1) X IO^^Mq. We plot, in Figure 0 the aver¬ 
age number of “building blocks” versus lookback time 
for our 58 Milky Way sized dark matter halos. To 
qualify as a “building block”, the merging sub-halo 
must be > IO^^Mq, add at least 4% to the mass of 
the halo, and be merged to the largest (central) halo 
by the next timestep. Thus, this plot summarises the 
history the build up of Milky Way sized halos (see 
for details). A peak merging period around 10 Gyrs 
ago, before a relatively quiescent final 7-8 Gyrs, is in¬ 
dicated by our plot, and provide support to the disk 
galaxy evolution scenario presented here. 


4. DISCUSSION 

The properties of disk galaxies have proven to be 
difficult to accommodate within the favoured Gold 
Dark Matter cosmological paradigm. Gomplex gas 
processes involving energy feedback, for example from 
quasars, super novae and stellar winds, and multiple 
gas phases, have proven difficult to model, yet these 
processes appear to be crucial in the formation and 
evolution of disk galaxies within the context of hier¬ 
archical structure assembly. Here we report on a sim¬ 
ulated disk galaxy, and match to recent observations 
of local and distant (z ^ 1) disk galaxies. In form¬ 
ing our simulated galaxy, feedback and gas physics 
proved crucial ingredients in forming a realistic disk 
galaxy. We now have have a suite of four simulated 
disk galaxies, and the galaxy presented here will form 
an extensive study of disk galaxy formation and evo¬ 
lution (i & 0 ) as we build up a statistical sample 
with which to interpret observations. 
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Figure 5: We zoom in to axis of 40 kpc, and show the /-band luminosity plots of our simulated galaxy, both edge-on 
(upper panels) and face on (lower panels). The flattened nature of the galaxy is apparent by the panel. The galaxy 
evolves relatively quiescently during the final 8 Gyrs. 
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Figure 6: /-band luminosity of the simulated galaxy 
versus radius, in the galactic plane, taken at three 
redshifts. The earliest timestep {z = 0.72, red) is taken 
after significant mergers have finished. An intermediate 
timestep (z = 0.46, yellow) and present (blue) are also 
shown. The exponential nature of the disks between 
5 — 15 kpc allows scale-lengths to be derived. 


The simulated disk galaxy is shown to have a thick 
disk component. This is evidenced by the velocity 
dispersion versus age relation for solar neighbourhood 


Figure 7: /-band luminosity versus height above the 
plane. Symbols as for Figure |H| Exponentials are fitted 
between 0.5 and 1.5 kpc to derive scale-heights. 


stars (Figure^. Our simulations indicate that the 
thick disk is created in an epoch of of multiple merg¬ 
ers of gas rich buildiM blocks, during which a cen¬ 
tral galaxy is formed Q, and it was shown that such 
a scenario is consistent with present observations of 
the Galactic and extragalactic thick disks. Recently, 
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Figure 8: The average number of building blocks versus 
lookback time (see text for details). 


observations of four edge-on galaxies explored proper¬ 
ties of their stellar populations as a function of height 
0. All showed thick disks composed of red stel¬ 
lar populations, which are old and metal rich. V—I 
colour gradients were seen to be zero or slighlty pos¬ 
itive. These observations do not support thick disks 
which form through dissipative disk formation, accre¬ 
tion from shredded satllites, or slow heating of the 
thin disk. Heating through galaxy interactions more 
than 3 Gyrs ago are not ruled out. In Figures 01 &E1 
we show the metallicity and age of the thick disk stars 
in our simulated galaxy. We see that our thick disk 
formation scenario results in stellar populations which 
can match the observations. Although accretion from 
galactic “building blocks” are involved in our model, 
these accreted building blocks have high gas fraction. 
Astration during and after murders, as well as ongo¬ 
ing disspiative infall of gas, eliminates differences in 
metallicity of original and accreting material. Thus 
little metallicity gradient is seen in our simulated thick 
disk. The rapid star formation in this epoch ensures 
that thick disk stars are almost exclusively old, and a 
lack of gradient with height. 

In a study of 34 edge-on disk galaxies in the Hubble 
Deep Field 0, galaxies out to redshift^l were shown 
to have scale parameters similar to local disk galax¬ 
ies. These parameters are characterised by the ratio of 
scale-height to scale-length {hz/h\). Disks at redshift 
z 1 are smaller, in absolute value, than present day 
disk galaxies. However, their flatness ratio, h^/hx, is 
only slightly larger, by a factor ^1.5. We examine our 
simulated disk galaxy from a time after the thick disk 


formation epoch, when a late type galaxy has clearly 
formed. We derive a value of h^/hx which is larger by a 
factor ^ 1.2 in our simulated galaxy at 2 ; = 0.72 com¬ 
pared with the z = t). Thus, even though a chaotic 
period of merging which led to the emrergence of the 
thick disk was involved in the birth of our simulated 
disk galaxy, once a disk emerges, the global structure 
does not evolve to a large extent. We found the disk 
does get larger and slighlty flatter. 
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